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Using the axially-symmetric time-dependent Gross-Pitaevskii equation we study the Josephson 
oscillation of an attractive Bose-Einstein condensate (BEC) in a one-dimensional periodic optical- 
lattice potential. We find that the Josephson frequency is virtually independent of the number 
of atoms in the BEC and of the inter-atomic interaction (attractive or repulsive). We study the 
dependence of Josephson frequency on the laser wave length and the strength of the optical-lattice 
potential. For a fixed laser wave length (795 nm), the Josephson frequency decreases with increasing 
strength as found in the experiment of Cataliotti et al. [Science 293, 843 (2001)]. For a fixed 
strength, the Josephson frequency remains essentially unchanged for a reasonable variation of laser 
wave length around 800 nm. However, the Josephson oscillation is disrupted with the increase of 
laser wave length beyond 2000 nm leading to a collapse of a sufficiently attractive BEC. These 
features of Josephson oscillation can be tested experimentally with present set ups. 



PACS numbers: 03.75.Lm,03.75.Kk 



I. INTRODUCTION 

The observation of an oscillating Josephson current 
across an one-dimensional array of potential wells, usu- 
ally generated by a standing-wave laser field and com- 
monly known as an optical-lattice potential, in a trapped 
cigar-shaped Bose-Einstein condensate (BEC) by Catal- 
iotti et al. 1\ was the first manifestation of this phe- 
nomenon in trapped neutral bosons. Until then the 
Josephson effect was confirmed in superconductors with 
charged electrons and in liquid helium |2(. 

In the experiment of Cataliotti et al. 1] a Joseph- 
son oscillation was initiated in a repulsive 87 Rb BEC 
formed in a one-dimensional optical-lattice plus axially- 
symmetric harmonic traps by suddenly displacing the 
harmonic trap along the axial direction. They investi- 
gated the evolution of the Josephson frequency with the 
strength of the optical-lattice potential and found that, 
for a fixed laser wave length, the Josephson frequency 
reduces with the increase of the strength. 

There have been previous theoretical studies of Joseph- 
son oscillation to explain different features of the experi- 
ment of Cataliotti et al. [l| using the numerical solution 
of the time-dependent mean- field Gross-Pitaevskii (GP) 
equation 0. One 0,0,13 an d three-dimensional mod- 
els based on the GP equation have previously been used 
to study different aspects of Josephson oscillation. There 
have been other theoretical studies of Josephson oscilla- 
tion in a trapped BEC 0,HH using different approaches 
and under different conditions. 

In this paper we address several interesting features of 
Josephson oscillation in a trapped BEC as in the exper- 
iment of Cataliotti et al. 0]. How does the Josephson 
frequency evolve with the number of atoms and atomic 
interaction (attractive or repulsive) of the condensate? 
How does the Josephson frequency depend on the in- 
tensity (strength) and wave length of the standing-wave 
laser beam used to create the periodic optical-lattice po- 
tential? What are the limits on the optical-lattice poten- 



tial parameters for allowing a Josephson oscillation? Is it 
possible to have a collapse in an attractive BEC induced 
by a Josephson oscillation? 

In the actual experiment and a previous three- 
dimensional simulation [|| the extraction of Josephson 
frequency was rather indirect. The BEC in the optical- 
lattice potential was allowed to expand and form an inter- 
ference pattern of two bright pieces moving in opposite 
directions from a central piece after a certain interval 
of time. The position of the original BEC executing a 
Josephson oscillation in the periodic optical-lattice po- 
tential was predicted from the position of the observed 
interference pattern. In this indirect fashion the Joseph- 
son oscillation was identified and its frequency extracted. 
Although, because of the very small size of the BEC, this 
seems to be the only means for experimental observation, 
we employ a direct approach for the study of Josephson 
oscillation in the present numerical simulation, which is 
more accurate and easier to implement than the indirect 
approach employed in Refs. 0,13- In present simulation 
we directly follow the oscillating condensate without re- 
sorting to an expansion. 

There have been experimental ^(J ancl theoretical 
studies of the collapse of a 85 Rb BEC after suddenly turn- 
ing the atomic interaction from repulsive to attractive by 
manipulating an external magnetic field near a Feshbach 
resonance [l2|| . There has also been a study of induced 
collapse of an attractive BEC trapped by a system of 
laser beams in the presence of small fluctuations of the 
laser intensity '13] . In this paper, to the best of our 
knowledge, we study for the first time the possibility of 
collapse of an attractive BEC induced by Josephson os- 
cillation. In view of the routine experiments |14L lla ] on 
BECs formed in optical-lattice potentials, it should be 
possible to verify the predictions of the present study on 
the Josephson oscillation of an attractive BEC. 

In Sec. II we present the mean-field model based on the 
axially-symmetric time-dependent nonlinear GP equa- 
tion. In Sec. Ill we present the numerical results and 
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finally, in Sec. IV we present a summary of our study. 



II. MEAN-FIELD MODEL 

The time-dependent BEC wave function ^(r;t) at po- 
sition r and time t is described by the following mean- 
field nonlinear GP equation jjj 



ti'V 2 
2m 



+ V(r)+gN\*(r;t)f 



*(r; t)=0,(l) 



where m is the mass and N the number of atoms in the 
condensate, g — AitTi 2 a/m the strength of inter-atomic 
interaction, with a the atomic scattering length. In the 
presence of the combined axially-symmetric and periodic 
optical-lattice potentials V(r) = \muj 2 {p 2 + v 2 y 2 ) + V ^i 
where u> is the angular frequency of the harmonic poten- 
tial in the radial direction p, vlo that in the axial direction 
z, with v the aspect ratio, and = s E r cos 2 (k^z) is 
the optical-lattice potential created with the standing- 
wave laser field of wavelength A = 795 nm, as in the ex- 
periment of Ref. |H, with E R = h 2 k 2 L /(2m), k L = 2-k/X, 
and s (< 12) the strength. The normalization condition 
is j dr\V{r;t)\ 2 = 1. 

In the axially-symmetric configuration, the wave func- 
tion can be written as \t(r, t) — ip(p,z,t), where < 
p < oo is the radial variable and — oo < z < oo is the 
axial variable. Now transforming to dimcnsionlcss vari- 
ables p — \/2p/l, z — \/2z/l, t — tui, I = yjh/{muj) 1 and 



ip(p,z;r) = pyl 3 /V8i/j(p, z;t), Eq. JIJ becomes [If 
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where the nonlinearity parameter n — Na/l. In 
terms of the one-dimensional probability P(z,t) = 2tt- 
J dp\(p(p, z, t) | 2 / p, the normalization of the wave func- 
tion is given by dzP(z, t) = 1. 

The experiment of Cataliotti et al. was performed 
with repulsive 87 Rb atoms in the hyperfine state F — 
1,ttif = — 1 • The axial and radial trap frequencies were 
vuj = 2ir x 9 Hz and uo = 2tt x 92 Hz, respectively, with 
v = 9/92. In the present simulation we consider attrac- 
tive 85 Rb atoms in the hyperfine state F = 2, rriF = —2 
with mass quite close to 87 Rb atoms. The triplet scat- 
tering length ar for the 85 Rb- 85 Rb system is negative 
corresponding to attraction: — 500ao < ax < — 300ao 
with a = 0.5292 A the Bohr radius 0. Also, in case 
of 85 Rb the atomic interaction can be modified using a 
Feshbach resonance 01 ■ So the 85 Rb- 85 Rb system is the 
most suitable for studying the dependence of Josephson 
oscillation on atomic interaction. 



For 85 R b the harmonic-oscillator length / = 
yhj (mui) — 1.14 yitm and the present dimensionless 
length unit corresponds to Z/V2 = 0.8 pm. The present 
dimensionless time unit corresponds to ui^ 1 — l/(27rx92) 
s = 1.73 ms. Although we perform the calculation in di- 
mensionless units using Eq. (J2J), we present the results 
in actual physical units using these conversion factors 
consistent with 85 Rb atoms. In terms of the dimension- 
less laser wave length Ao = y2X/l ~ 1, the dimension- 
less standing-wave energy parameter Er/ (Tluj) = Att 2 /Xq. 
Hence in dimensionless unit the periodic optical-lattice 
potential of Eq. J5J is 



^opt E R 2 

— — = s— cos (k L z) 
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where we have introduced a parameter a to manipulate 
the laser wave length. Normally, a — 1 corresponds to 
the experimental situation: A = 795 nm. A general value 
for a simulates a laser wave length of aX. 

We solve Eq. J3J) numerically using a split-step time- 
iteration method with the Crank-Nicholson discretization 
scheme described recently 0|. We discretize the GP 
equation typically with time step 0.001 and space step 0.1 
spanning p from to 7 pm and z from —70 /im to 70 /im, 
although, sometimes we used smaller steps for obtain- 
ing convergence. The time iteration is started with the 
known harmonic oscillator solution of Eq. J2J for n = 0: 
<p(p,z) = [^/(Stt 3 )] 1 / 4 p e -<,P 2 +v* 2 )/i with chemical po- 
tential p, = (1 + v/2) [TaL For a typical cigar-shaped 
condensate with ^~0.1 l| p ~ 1 is much smaller than 
the typical depth of the optical-lattice potential wells 
E r /(Huj) = 47r 2 /Ao ~ 40 so that p « Er/(Hu>) and the 
passage of condensate atoms from one well to other can 
only proceed through quantum tunneling. The nonlin- 
earity n as well as the optical-lattice potential parameter 
s are slowly increased by equal amounts in lOOOOn steps 
of time iteration until the desired value of nonlinearity 
and optical-lattice potentials are attained. Then, with- 
out changing any parameter, the solution so obtained is 
iterated 50 000 times so that a stable solution is obtained 
independent of the initial input and time and space steps. 
The solution then corresponds to the bound BEC under 
the joint action of the harmonic and periodic optical- 
lattice potentials. 



III. NUMERICAL RESULTS 

First we consider an attractive BEC formed in the com- 
bined harmonic and periodic optical-lattice potentials for 
a specific attractive (negative) nonlinearity n. Such a 
BEC is stable for \n\ less than a critical value depending 
on the parameters of the optical-lattice potential . We 
study the formation of a BEC in the combined harmonic 
and optical-lattice potentials for a = 1 in Eq. © for a 
range of values of s. Depending on the value of a, the 
BEC could be unstable for \n\ > 0.4 [n| and for a = 1 
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FIG. 1: The profile of the BEC wave function \ip(p, z) \ formed 
in the joint optical-lattice and harmonic potentials for n — 
—0.3, s = 4, and a = 1. In the off-set the one- dimensional 
probability P(z) of this wave function vs. z is plotted for 4 /an 
> z > — 4 pm which clearly exhibits the maxima and minima 
of the condensate wave function along the axial direction. 



we consider n = —0.3. However, for a larger a, the BEC 
is stable for a larger value of \n\ 0] (see, for example, 
the stable condensate studied in Fig. 5 for n — —0.5 in 
the following for a = 3). In Fig. 1 the plot of \ip(p, z)\ 
vs. p and z illustrates the profile of the BEC bound state 
for n = —0.3 and s = 4. The BEC is in the form of nar- 
row slices with the optical-lattice barriers separating the 
slices. The large number of maxima and minima due to 
the optical-lattice potential is not clearly visible in this 
plot. The maxima and minima in the axial direction are 
clear in the plot of the one-dimensional probability P(z) 
vs. z in the off-set of Fig. 1 for 4 /im > z > —4 /im. In 
this interval of z, there are 16 wells of the optical-lattice 
potential and as many maxima and minima in P(y,t). 

Next we consider an oscillating BEC in the combined 
harmonic and periodic optical-lattice potentials. If we 
suddenly displace the harmonic trap along the lattice 
axis by a small distance after the formation of the BEC 
in the combined potentials, the condensate will acquire 
a potential energy, be out of equilibrium and start to os- 
cillate. As the height of the potential-well barriers of the 
optical-lattice potential is much larger than the energy 
of the system, the atoms in the condensate will move by 
tunneling through the potential barriers. This fluctuat- 
ing transfer of atoms across the potential barriers is due 
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FIG. 2: (Color online) Contour plot of the one-dimensional 
probability density P(z, t) of the attractive BEC with a = 1 
and n = —0.3 executing a Josephson oscillation when the 
harmonic trap is suddenly displaced along the axial direction 
through a distance 16 pm at t = for optical-lattice strength 
s = (a) 1, (b) 2, (c) 3, (d) 4, (e) 5 and (f) 6. The period of 
Josephson oscillation can be obtained from these plots. 



to Josephson effect in a neutral quantum liquid. 

With the attractive BEC of Fig. 1 we next study its 
Josephson oscillation when the harmonic potential is sud- 
denly displaced along the axial direction by 16 /im. The 
BEC now acquires an added potential energy which it 
can use to execute a Josephson oscillation along the ax- 
ial direction. The Josephson oscillation is best studied 
numerically from the contour plot of the one-dimensional 
probability P(z,t) vs z and t exhibited in Figs. 2. These 
plots clearly show the central position of the BEC along 
the axial z direction. In the present simulation we take 
a = 1 and different values of the optical-lattice strength 
s. These contour plots are very useful to find the Joseph- 
son frequencies. From Figs. 2 the periods of Josephson 
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FIG. 3: (Color online) The frequency of the atomic current 
in the array of Josephson junctions as a function of optical- 
lattice strength s: • with error bar — experiment of Cataliotti 
el al. p| with repulsive 87 Rb atoms; A — tight binding calcu- 
lation taken from Ref. 1]; * — simulation with repulsive 87 Rb 
atoms taken from Ref. |(j; full line - present calculation with 
attractive atoms; dashed line - hydro-dynamical calculation 
with repulsive atoms taken from Ref. 0. 



oscillation are easily read off and the frequencies calcu- 
lated for different optical-lattice strengths s. 

In Fig. 3 we plot the Josephson frequencies vs. optical- 
lattice strength s from different studies for a — I. Specif- 
ically, in addition to the present calculation with at- 
tractive atoms we also show (1) the experimental fre- 
quencies and corresponding tight-binding calculations for 
repulsive atoms of Cataliotti et al. 1], (2) the three- 
dimensional simulation in the repulsive case from Ref. 
0, and (3) the hydro-dynamical calculation with repul- 
sive atoms from Ref. 0. 

The most interesting conclusion from Fig. 3 is that the 
present frequencies of three-dimensional simulation for an 
attractive BEC are practically the same as the frequen- 
cies of the repulsive BEC of Ref. Q • Hence the Joseph- 
son frequencies are independent of, or at best weakly de- 
pendent on, the number of atoms in the condensate and 
nature of interaction (repulsive or attractive). This is 
because in Eq. J2J) the information about the number of 
atoms and atomic interaction is contained solely in the 
nonlinearity parameter n = Na/l. For s = the optical- 
lattice potential is absent and the condensate executes 
free oscillation with the frequency of the axial potential 
(9 Hz) and the three-dimensional calculations by this au- 
thor converges to this value in the s — > limit for both 
attractive and repulsive interactions. 

However, there are some discrepancies between the 
present three-dimensional calculation and the tight- 
binding calculation of Ref. as well as between the 
present calculation and the hydro-dynamical calculation 



of Ref. 7] for small values of optical-lattice strength s. 
The Josephson frequencies of Refs. 0, Q do not seem 
to lead to the same result (9 Hz) as the present calcula- 
tion in the low s limit. The reason for this could be the 
possible use of a slightly different axial trap frequency in 
these studies. For large values of optical-lattice strength 
s, there are no experimental results. In this domain, the 
present frequencies ^gree well with the hydro-dynamical 
calculation of Ref. [?J ■ But both these theoretical results 
differ from the tight-binding calculation of Ref. . The 
reason for this could be the use of a variational Gaussian 
ansatz in Ref. 0. 

Now we study the effect of the variation of laser wave 
length on the Josephson frequency for a fixed strength 
s of the optical-lattice potential. To implement it we 
vary the parameter a of the optical-lattice potential 
which effectively varies the wave length of the laser. The 
Josephson frequency is found to be weakly dependent on 
laser wave length for 1600 nm > A > 600 nm leading to 
frequencies 7.1 ± 0.3 Hz. If Josephson frequency were 
solely determined by the effective strength of the optical- 
lattice potential (47t 2 /Aq) x (s/a 2 ) of Eq. (3), then from 
Fig. 3 a much larger variation of frequency would be 
expected for a variation of A in the range 600 nm to 1600 
nm. The weak dependence of Josephson frequency on a 
suggests that in addition to the strength of the optical- 
lattice potential, the cosine term in Eq. (3) also plays a 
decisive role in determining the Josephson frequency. 

Next, we consider the effect of increasing the laser wave 
length beyond A = 2000 nm. The condensate in the 
periodic optical-lattice potential is in the form of nar- 
row slices with the optical-lattice barriers separating the 
slices (see Fig. 1 and related description) . These slices of 
the BEC have to tunnel through the optical-lattice bar- 
riers during Josephson oscillation. The increase in laser 
wavelength corresponds to an increase in the widths of 
the slices of the BEC and optical-lattice barriers and to 
execute a Josephson oscillation, larger pieces of conden- 
sates have to tunnel freely through wider optical-lattice 
barriers. 

With the increase of the laser wave length beyond a 
critical value, the larger pieces of condensates cannot 
freely tunnel through wider barriers and there is actually 
a disruption of Josephson oscillation. This is illustrated 
in Figs. 4 through a contour plot of one-dimensional 
probability density P(z,t) with a = 3 and 4 when the 
harmonic potential is suddenly displaced along the axial 
direction through 16 /xm at t = for s ~ 4 for a small 
nonlinearity |n| = 0.1. The clean Josephson oscillation 
tracks seen in Figs. 2 now disappear with the increase 
of laser wave length as can be seen in Figs. 4 (a) and 
(b) where we show the contour of the one-dimensional 
probability density P(z, t) for a = 3 and 4. From Figs. 
2 (d), 4 (a), and 4 (b) we see that for s — 4 there is a 
clean Josephson oscillation track for a = 1. The same 
is true for a = 2. However, with an increase of a to 3 
the trail left by the condensate is not so clean showing 
the beginning of breakdown of Josephson oscillation in 
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FIG. 4: (Color online) Contour plot of the one-dimensional 
probability density P(z,t) of the BEC with s = 4 and n = 
—0.1 when the harmonic potential is suddenly displaced along 
the axial direction through a distance 16 y,m at t = for laser 
wave length parameter a = (a) 3 and (b) 4. The disruption 
of Josephson oscillation has already started for a — 3 and is 
complete for a = 4. 



Fig. 4 (a). There is a complete breakdown (absence) of 
Josephson oscillation in Fig. 4 (b) for a = 4. 

There has been a different account of breakdown of 
Josephson oscillation in a neutral quantum fluid in Ref. 
1 2 ("| when the harmonic potential was displaced by a very 
large distance (> 130 /im) along the periodic optical- 
lattice potential. It was suggested that the breakdown 



of Josephson oscillation in that study was due to a loss 
of phase coherence in the condensate due to a classical 
transition from a superfluid to an "insulator" resulting 
in a modulational instability. Other mechanisms for the 
loss of phase coherence in a condensate have also been 
studied |2l|. The loss of phase coherence considered in 
all these investigations [20l l2l| originated from a clas- 
sical superfluid to insulator transition, different from a 
quantum transition of a superfluid to a Mott insulator 
observed in Ref. [l4j . 

From Fig. 3 we find that for a fixed laser wave length, 
the Josephson frequency reduces with increasing optical- 
lattice strength, which seems plausible. From Figs. 4 
we find that for a fixed optical-lattice strength with an 
increase in laser wave length the Josephson oscillation is 
strongly attenuated. As the Josephson oscillation takes 
place by quantum tunneling in atomic superfluid through 
the optical-lattice barriers, this oscillation is bound to be 
reduced as the height and width of the optical-lattice 
barriers are increased. This reduction of Josephson oscil- 
lation with the increase of the parameters s or a will first 
manifest in a reduction of Josephson frequency and even- 
tually to a disappearance of Josephson oscillation. This 
reduction in frequency is explicit in Fig. 3. A breakdown 
of Josephson oscillation with increasing a is found in Fig. 
4. The possible breakdown of Josephson oscillation with 
large s for a = 1 is not investigated in this paper. 

Finally, we investigate if the above disruption of 
Josephson oscillation for a large a could provoke a col- 
lapse in an attractive BEC. We find in our numerical 
study that in the cases where there is a clean Josephson 
oscillation there is no collapse in an attractive conden- 
sate. However, for large a, when there is a disruption of 
Josephson oscillation, e.g., as in Figs. 4, one could have 
a collapse in an attractive condensate with sufficiently 
strong nonlincarity. 

To demonstrate the above claim we performed the 
simulation of Josephson oscillation using the parameters 
used in Figs. 2 for attractive nonlinearity close to the 
critical nonlinearity for collapse 0|. In the actual con- 
ditions of Figs. 2, the Josephson oscillation is always 
perfect and there is no collapse in the condensate pro- 
voked by this oscillation. However, the scenario changes 
as the parameter a is increased for a large nonlincar- 
ity \n\. To illustrate this we performed calculation with 
optical-lattice strength s = 4, a — 3, and nonlinearity 
n = —0.5. The stationary BEC wave function in this case 
is localized and stable as shown in Fig. 5 (a) at t — 0. A 
Josephson oscillation is then initiated by displacing the 
harmonic potential along the axial direction through a 
distance of 16 /an. Soon after this displacement, due to 
the onset of collapse there is no stationary solution to 
the GP equation. Numerically, the nonexistence of a sta- 
tionary solution to the GP equation and the subsequent 
nonequilibrium dynamics of the condensate manifests in 
a destruction (delocalization) of the localized stationary 
wave function at t = and a generation of a dynami- 
cal nonequilibrium wave function spread over the entire 
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The investigations of Refs. 0, @ also provides an- 
swer to some of the questions raised in the present pa- 
per, e.g., how does the Josephson oscillation depend on 
the nonlinearity parameter n and the optical-lattice po- 
tential strength s. Using a hydro-dynamical model, the 
authors of Ref. show for an elongated or cigar-shaped 
confinement (y « 1) that the frequency of oscillation 
lu z = vuj in a pure harmonic potential gets changed in 
the presence of the optical-lattice potential according to 
^opt = UJz \/ m l m * i where m is the atomic mass and m* 
is an effective mass. A good estimate of m/m* can be 
made by neglecting both the magnetic trapping and the 
atom-atom interaction from the spectrum of the linear 
Schrodinger equation for the one-dimensional periodic 
optical-lattice potential. In this approximation m/m* 
is a universal function of the parameter s, calculated in 
Ref. 0. This dependence of m/m* on s has been used 
to calculate the dependence of Josephson frequency on s 
as quoted in Fig. 3 of this paper. The quantity m/m* 
also depends on nonlinearity n and harmonic oscillator 
frequency to. However, the dependence of m/m* on n 
and lu is weak for small values of these parameters 0. 
In Ref. the authors study the dependence of m/m* 
on nonlinearity n in a model and demonstrate the weak 
dependence for small and medium values of n. Conse- 
quently, the independence of the Josephson frequency on 
the nonlinearity parameter n as well as the dependence 
of this frequency on the parameter s as noted in Fig. 3 
seems to be consistent with Refs. 171 However, it 
should be noted that the studies of Refs. did not 
consider Josephson oscillation for attractive interaction 
and the possibility of collapse as in this paper. Hence 
the conclusions of Refs. 0,13 are only valid for repulsive 
atomic interaction and the present study extends them to 
the case of attractive atomic interactions. Moreover, the 
present study is based on the three-dimensional mean- 
field microscopic GP equations, whereas those of Refs. 
Q are based on hydro-dynamical and one-dimensional 
mean-field models of the condensate. 



FIG. 5: (a) The profile of the wave function \ip(p, y)\ formed 
in the joint optical-lattice and harmonic potentials for s = 4, 
a = 3, and n = —0.5 at t = 0. (b) The profile of the wave 
function of (a) at t = 16 /im after an axial displacement of 
the harmonic potential through 16 /im at t — 0. 



spatial domain of discretization as seen in Fig. 5 (b) at 
t = 16 ms. 

When there is a free Josephson tunneling, the conden- 
sate oscillates as a whole without significant modification 
of shape or local density. However, when there is a dis- 
ruption of Josephson tunneling and the BEC is given 
additional potential energy to oscillate by displacing the 
harmonic potential, the condensate is locally squeezed to 
a smaller region in space resulting in a very large den- 
sity and thus inducing a collapse of the condensate as 
illustrated in Figs. 5. 



IV. SUMMARY AND CONCLUSION 

We have performed numerical simulation based on an 
axially-symmctric mean-field GP equation to address dif- 
ferent aspects of Josephson oscillation of an attractive 
BEC in a combined harmonic and periodic optical-lattice 
potentials. The latter potential is created by a standing- 
wave laser beam. We find that the Josephson frequency 
is virtually independent of the number of atoms in the 
condensate and the nature of atomic interaction: attrac- 
tive or repulsive. The mean-field results for Josephson 
frequencies for different optical-lattice strengths are in 
qualitative agreement with experiment 0] and other cal- 
culations IE 0- However, there is some quantita- 
tive disagreement between different results. The Joseph- 
son oscillation continues with essentially unchanged fre- 
quency (7.1 ± 0.3 Hz) for a variation of laser wave length 
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in the region 1600 nm > A > 600 nm. However, as A is 
increased past 2000 nm there is a disruption of Joseph- 
son oscillation. For an attractive BEC with nonlinearity 
stronger than a critical value this disruption of Joseph- 
son oscillation could induce a collapse. These features of 
Josephson oscillation could be tested experimentally and 
will provide a test for the mean-field model. 
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